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Abstract 

It is possible to discuss the propagation of an electronic current through 
certain layered nanostructures modeling them as a collection of random 
one- dimensional interfaces, through which a coherent signal can be trans- 
mitted or reflected while being scattered at each interface. We present a 
simple model in which a persistent random walk ( the "t-r" model in 1-D) 
is used as a representation of the propagation of a signal in a medium 
with such random interfaces. 

In this model all the possible paths through the system leading to 
transmission or reflection can be enumerated in an expansion in the num- 
ber of loops described by the path . This expansion allows us to conduct a 
statistical analysis of the length of the paths for different geometries and 
boundary conditions and understand their scaling with the size of the sys- 
tem. By tuning the parameters of the model it is possible to interpolate 
smoothly between the ballistic and the diffusive regimes of propagation. 
An extension of this model to higher dimensions is presented. We show 
Monte Carlo simulations that support the theoretical results obtained. 



1 Introduction 

The seminal work of Anderson||] raising the possibility that disorder can lead to 
non-diffusive behavior (the so called localized regime) refocused the attention of 
the Physics community on the problem of the propagation of waves in disordered 
systems. In the last two decades new theoretical ideas (like the scaling the- 
ory of localization|Q, weak localization!^, universal conductance fluctuations Q 
and Wigner dwelling times ||^) were advanced, and a new field (soon called 
Mesoscopic Physics) emerged. It reached and influenced many experimental 



areas, among them electronic systems[||, microwaves|0, optics|^, acoustics]^, 
geophysics pO|, laser physics [|ll| , medical physics Q and atomic physics Q. 

It has become an extremely important problem in this field to understand 
what should be the signature of the propagation of a signal in the different 
regimes (ballistic, diffusive, localized) since concomitant phenomena, like ab- 
sorption can complicate the interpretation of experimental results . It is for 
that reason that theoretical analyses of the characteristics of the propagation, 
and in particular its statistical properties are of great interests, since 
those properties have become recently experimentally accessible0[^. 

When the inelastic scattering length in a system is large compared to its 
size, the wave propagates coherently in the sense of its phase being preserved 
while its direction is randomized by elastic scattering processes with the impu- 
rities constituting the random medium. It is natural, then, to establish con- 
nections between the coherent propagation of a wave and the statistical theory 
of random walks ||l8|, where the continuum limit is known to lead to diffusive 
theory [|9j. This connection is relevant to the propagation of an electronic 
current through certain layered nanostructures, since they can be modeled as 
a collection of random one-dimensional interfaces, through which a coherent 
signal can be transmitted or reflected while being scattered at each interface. 

The connection between the effect of the disorder of the propagating medium 
and the statistical randomness of the mathematical description of the diffusion 
has not been satisfactorily clarified yet. One aspect that did not receive enough 
attention in this approach is the fact that, since the scattering with the impu- 
rities is not isotropic and the cross section is normally enhanced in the forward 
direction, the statistical jumps are not independent from each other. In other 
words the probability distribution of each step of these persistent random walks 
(PRW) is dependent on the previous step|^. The PRW are more difficult 
to study than the standard random walks, where the probability distribution 
of each step is independent of what happened in the previous steps. This can 
be seen from the fact that, for example, while a standard random walk can 
be easily mapped into a chain of non interacting spins (in the presence of a 
magnetic field, if the random walk is biased), the same mapping for the PRW 
leads to a chain of nearest-neighbor interacting spins (the full Ising model for 
the case of a 1-D PRW) . In particular in this work we will explore the 
connection between the propagation of a wave and the 1-D version of the PRW 
called the "t-r" model in which a particle moves with probability t in the 
same direction as in the previous step and with probability r reverses direction 
(t + r = 1). This type of model is more suited than the standard random walk 
model to explore the quasiballistic regime, which becomes more important in 
transport phenomena in nanostructures. 



2 The model 

We will consider a model of the propagation of a wave first introduced for nu- 
merical purposes by Edrei, Kaveh and Shapirop2|, and later applied extensively 
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by Vanneste et aZ[Q and SebbahQ. Each site is represented by a scattering 
matrix S of dimension 2D x 2D [D being the dimensionality) connecting the 
2D incoming amplitudes at time t, with the 2D outgoing ones at time t -\- 1. 
These amplitudes can be regarded as residing in the bonds (2 joining each site 
to its nearest neighbor, one in each direction). These bonds represent free prop- 
agation between the sites, and their phase can be included in the S matrix. The 
wave function at time t is, thus, determined by the complex numbers associ- 
ated with each of the bonds at that time. Numerically this has the advantage 
of letting the stationary wave being built step by step in time, from an initial 
input without solving a huge diagonalization problem. The matrices S repre- 
senting each scatterer are unitary and symmetric reflecting, respectively, the 
energy conservation and time reversal symmetry inherent to the problem. If we 
adopt the simplifying assumption that the scatterers themselves are symmetric, 
these matrices simplify even further, since the following relations can be easily 
obtained for their elements 

^-t^ + 2{D = 1 (1) 
rt cos[Lpr - (fit] + {D - = 
rd cos[(pr ~ (fid] +tdcos[ipt — (pd] + {D — 2)d^ = 

where {r, </9r}, {t, ipt} and {d, Lpd} are the amplitudes and phases of the elements 
for reflection, transmission, and turning (in ID, d = 0). The matrices S have, 
thus, only two (in ID) or three (in D>1) independent parameters. For example, 
in ID the scattering matrix representing one scatterer can be parameterized as: 

exp[/(^] {r^+t'^) = l (2) 

In this simple version of the model the distances are measured in units of the 
mean free path (the distance between scattering elements) and the phase velocity 
of the wave is one. The effect of the disorder in this model has been considered 
both analytically p5| and numerically p6| , by assuming a distribution of values 
for the variable r, linked to the reflectivity of the scatterers. Alternatively, one 
can consider the effect of disorder in the variable ip, more closely linked to the 
distance between scatterers [p7[ . 

We will be interested in considering a signal incident from the left on a 
system composed of N such scatterers, and in monitoring what is reflected or 
transmitted through it at different times, i.e. the output to the left of scatterer 
1 and to the right of scatterer N, respectively. It is important to observe that 
while calculating the response of the system to an applied pulse implies adding 
all the paths of a certain length (that will interfere by arriving simultaneously to 
the boundary) , the problem of determining the response to a continuous wave is 
simpler, since in the stationary state it is sufhcient to add all the transmitting 
or reflecting paths as long as the phase factors associated with the propagation 
through each scatterer, exp[/fc], are included, k represents here the wavevector 
of the incident wave. 
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Figure 1: (a) Three transmitting paths through a system of two scatterers, 
showing no loops, one loop of length one, and two loops of length one, re- 
spectively, (b) Three transmitting paths through a system of three scatterers, 
showing no loops, one loop of length two, and the same loop of length two 
"dressed" by two loops of length one, respectively. 



3 Loop expansion 

The amplitude factor for any particular path reaching the boundaries of the 
system can be easily written, and we can sum all the possible paths by perform- 
ing a loop expansion. For example, for a system of two scatterers the ampli- 
tudes associated with the simple paths illustrated in figure 1(a) are, respectively, 
(/fexp[/A:])2 , (/texp[/fc])2(rexp[/A:])2 and {I te^p[I k]f{r eiip[I k])^; the to- 
tal amplitude for transmission through the two scatterers, when all possible 
paths are considered, becomes then: 

(/texp[Jfc])^ _ (/texp[/fc])2 
' l-(rexp[/A:])2 ~ I - h ^ ' 

For the three scatterers in figure 1(b), the amplitude factors associated with 
the paths shown are {I t exp[I k])^ , (I tcxp[I k])^{r exp[I k])'^{I texp[I k])"^ and 
{Itcxj>[Ik]f{rcxp[Ik]f{Itcxp[Ik])^ (rexp[/fc])4, respectively The "bare" 
loop of length 2 has an amplitude factor of I2 = (r cxp[I k])'^ {I t cxp[I k])'^ , and 
when fully "dressed" by all the possible loops of length 1 becomes: 

~ ^ (rcxp[/ k])-[Ilcxi,[I k]f ^ 
' (l-(rexp[7fc])2)2 ~(i_/i)2 

We see that the factor associated with a "bare" loop of length n is simply 
In = (rexp[7 A;])2(7texp[/ A;])2("~i). In terms of it we can obtain a recursive 
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expression for the general "dressed" loop of order n 



li — li ; I2 — — — ; ^3 = — — 



/„ = ^ ^ (4) 

(l_|,)2(n-l)(l_^2)2(n-2) 

It is possible now to add all the paths through a chain of scatterers, the 
amplitude of the transmission through the system, and express it in terms of 
the loops: 



(/f exp[//cl)2 (/iexp[/fc])3 
ii =7i exp[7fc] ; t2 = L ; = 



(1-/1) 

{It exp[7fc])" 



(l-Zi)("-i)(l-Z2)("-2) 



(5) 



As expected, In is proportional to tn^. A very similar expansion yields, for the 
reflection amplitudes: 

(/iexp[/fc])2(rexp[7A;]) 
ri=rexp[7fc] ; r2 = n + — — ^^-4 — — - 

(l-^i) 

(7texp[7fc])4 (rexp[7fc]) 
(1-/03(1-/;) 



r-3 = r2 + 



(7iexp[7A:])2("-i)(rexp[/fc]) 

I'n = 1~n-l + — — — — (6) 

(1 - /l)(2"-3)(l - Z2)(2n-5) (1 _ Z„_2)3 _ ;^_^) 

The expression for the amplitude of reflection through the system admits a 
clear interpretation in terms of a sum of all the paths that penetrate up to a 
certain depth in the system: ri corresponds to the paths that do not enter into 
the system at all, r2 — ri represents the paths that penetrate up to a depth 1, 
rz — r2 are the paths that penetrate up to a depth 2, etc. 

The connection of this expansion for the transmission and reflection ampli- 
tudes of the propagating wave with the PRW is provided by the fact that we 
can write exactly the same expansion in loops in the "t-r" model, for all the 
possible paths of a random walker leading to transmission or reflection. Here 
transmission means the random walker arriving to a point further to the right 
of position A'^, and reflection the walker reaching a point to the left of position 1. 
In particular, if we ignore the phase factors exp[7 k] and the imaginary unit in 
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front of t, the amplitude factors discussed above become the probabihties of the 
corresponding PRW paths, and Equations (^,(||) and (^) describe the "dressed" 
loops, probability of transmission and probability of reflection in the case of a 
PRW as well, since for a random walk the total probability of transmission 
(reflection) is simply the sum of the probabilities of each of the paths leading 
to transmission (reflection). Of course in each case, for a PRW the connection 
between t and r is not provided by (||), but by the simple expression t + r = l. 

It is easy to show now that the probability of transmission for a PRW be- 
comes, then: 



*"~l + (iV-l)r 

and the probability of reflection becomes simply r„ = 1 — i„ . Identifying N 
with the length of the system L, and taking into account that in this model 
< cos{-d) >= t — r, we can see that in the limit of large N the transmission 
tends to the diffusive limit t{L) — {I* /L). 

In the case of a wave, the transmission or reflection probabilities are given 
by the modulus square of the expressions calculated in (^ and (^ . This implies 
sums of the type: 

2 

Y^A, =^|A,r + 5]A,A* (8) 

i i i^j 

It is easy to see that, when adding the transmitting paths' amplitudes, if we 
consider the contribution of only the first term in the expression above (neglect- 
ing the correlations between the paths) we will reobtain (|^), where the roles 
of t and r are played by the transmission coefficient = T, and the reflection 
coefficient = R. It is precisely the overlap between paths (the second term) 
that leads to corrections to the diffusion theory results. 



4 Statistics of the paths 

In the case of the PRW, knowing the probabilities pi and the lengths U of all the 
paths both in transmission and reflection, we are now in a position to study some 
of their statistical properties. Defining a generating function Ftr (s) = ^ pi s'* , 
where we sum over all the transmitted paths, and a similar definition for Fref (s) 
where the sum is over all reflected paths, we can calculate the average length of 
the paths as 



< Z > = lim ^ ^'^} , 
s^i F{s) ' 



obtaining for the average length: 



l + {N-2)r+l{N-2){N-l)r^ 
< ^ ^ [l + (iV-l)r][l-r] 
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Ar+l(2iV-l)(iV-l)r 
< ^ l + (iV-l)r 

In the diffusive limit discussed above, these expressions reproduce the well 
known results of diffusive theoryQ: < I >tr (L) (L^/l*). Analogously, 
< I >ref [L) (2/3)i, independent of the mean free path. The expressions 

and ( p^ ) are exact, within the model, and remain valid in the diffusive as 
well as the ballistic limit r — > 0, and can be used to interpolate between them. 
This can be of interest in any situation in which the ballistic paths play an 
important role as in medical imaging. 

In a similar way one can express in terms of F(s) the variance of the path 
length distribution, in transmission or reflection 



obtaining the expressions: 



N(N^ - 1) r2 [N + hN - 2){2N - 1) r + ^{N - 1){N - 2){N - 3) 
vartr = — 



(l-r)2[l + (iV-l)r]2 



(11) 



""^-^ = 3(l-r)2 [l + (iV-l)rP ^'2) 

P(r) is 1 + (2A^ - 3)r + ^{UN^ - SON + 19)r^ + -^{2N - 1){N -1){N- 2)r^ 
here. These exact expressions, which again allow us to extrapolate between 
the diffusive and ballistic regimes, in the diffusive limit tend to {L"^/ l*)'^ and 
{L^ /I*) respectively. In the case of transmission the variance scales with the size 
of the system as the square of the average length (non-self averaging property) 
while in the case of reflection it scales even faster than it, exhibiting the well 
known noisy character of reflection [p8|. 

Figures 2, 3 and 4 exhibit the results of numerical simulations of a persistent 
random walker on a system of size N — 100. The transmission and reflection 
probabilities, the average length of the transmitted and reflected paths and the 
variance for the path lengths in transmission and reflection are represented as 
a function of the transmission parameter t. In all cases the agreement with the 
expression and (|) to (|l^ is excellent. 

An interesting problem is how to deflne the probabilistic weights to perform 
a similar statistical analysis of the length of the paths in the case of a wave, 
since the amplitudes are complex numbers with a phase, and the different paths 
interfere with each other. It makes sense to use the "probabilities" (in the 
sense of the moduli of the complex amplitudes squared) as those probabilistic 
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Figure 2: Fraction of transmitted paths (•) and of reflected paths (o) for different 
values of the transmission parameter t in a numerical simulation of up to 210,000 
random walkers on a system of size N = 100. The full lines correspond to 
equation (|^) for iioo , the probability of transmission, and rioo the probability 
of reflection. 



Figure 3: Average length of (a) the transmitted paths and (b) the reflected 
paths for the same simulation. The full lines correspond to the equations 
and ( 0) for N = 100. 
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Figure 4: Standard deviation of the distribution of the length of (a) the trans- 
mitted paths and (b) the reflected paths for the same simulation. The full lines 
correspond to the square root of equation (11) and equation (p^ for N = 100. 



weights. As the discussion at the end of section g shows, those weights describe 
the diffusive behavior of the system. Adopting this ansatz we can define again 
the generating function as in the case of the PRW, and we can reobtain all the 
results in expressions ^ to (|l2|), again replacing t and r by the transmission 
coefficient T and the reflection coefficient R, respectively. 



5 Higher dimensions 

All this analysis and the expansion in terms of loops can be easily generalized to 
a system in a higher number of dimensions. To flx ideas let us consider, instead 
of a row of N scatterers, a strip of M rows of N scatterers each, where the wave 
(or the random walker) is incident from left. For this system, the amplitudes 
of transmission or reflection through one column become (instead of t and r) 
the M X M matrices Ti and Ri , whose element (Ti)y is the amplitude for 
a walker incident from the left on row i to exit to the right on row j, and an 
analogous definition for (Ri)ij. These amplitudes are obtained by adding all 
the possible paths that lead from the input i to the output j and are confined 
to one single column of the medium. It is, then, not difficult to prove that 
the expansion in loops (Q) remains valid as long as we replace in it t and r 
by these matrices, and all the products are understood as products of matrices. 
Similarly, the transmission T^v and the reflection Rjv through the strip of length 
N are obtained through the same replacement from the expressions (|^) and (^ 
respectively. 

If the signal impinging on the system from the left is represented by an 
M-component vector I, the output vector O simply becomes, for the case of 



9 



Figure 5: Transmission through a system of length N—S as a function of t, for 
the cases of a 1-D chain, a strip of width AI — 2 (reflecting and absorbing 
boundary conditions) and for a 2-D periodic system. The two-dimensional sys- 
tems correspond to the case d = r. The inset illustrates the "snaking" effect for 
t 0. 



transmission 



Since in this model the matrices Ti and Ri have exactly the same symmetries, 
they will have the same eigenvectors; those will also be the eignevectors of the 
complicated matrices Tn and R^v, since they are expressed (through the loop 
expansion) as products of Ti and Ri. The elements of the matrix Tjv involved 
in (l^can be evaluated in terms of those eigenvectors and the eigenvalues of Ti 
and Ri through standard algebraic techniques. The problem is, then, essentially 
reduced to solving the eigenvalue problem of an M x M matrix. 

To conclude we will illustrate the treatment of the two dimensional problem 
by considering the particularly simple case of periodic boundary conditions in 
the lateral sides of the system (rows 1 and M). We will calculate the total 
transmission through a strip of length N. Then, 



where the stand for the eigenvalues of the matrix Tn, and the E" for the 
corresponding eigenvectors of the matrices Ti and Ri. From the symmetry of 




(13) 




(14) 
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the problem readily follows that one of the eigenvectors is Eq = -^[l^ 1,1, 1) 

and all the others cancel when summed over j in by orthogonality, leading 
to Tat = (Sjl*)- The eigenvalue can be simply obtained by replacing 
in the expansion or its equivalent for a PRW the coefficients t and r by the 
eigenvalues A° ^ and AJ ^ (associated with Eq) of the matrices Ti and Ri. 
Those eigenvalues are simply for any row i, and a similar expression 

for Ri. For the case of a PRW[p9|, (where the walker has a probability t of 
moving forward , a probability r of moving back, and a probability d of turning 
to either side, with t + r + 2d = 1), these eigenvalues can be obtained writing all 
the possible paths starting from position i that lead to transmission (reflection) 
through any row of a 1-column system: 



j 

The expansion in this expression corresponds to all paths with no reflection in 
the vertical direction, one reflection in the vertical direction, etc.). This is simply 
X°j, = t + d. Analogously we obtain A'j' ^ = r + d . The total transmission for 
a system of length N becomes, then: 



2 



2d^ 



-N 



(t + d) 



+ iN -l){r + d) 



Figure 5 illustrates the transmission through a strip of length 8, for periodic, 
absorbing and reflecting boundary conditions, and for a 1-D chain of the same 
length, in the specific case r = d. The inset illustrates the effect of "snaking", 
the existence of a non-zero total transmission even in the case of a zero forward 
transmission t at the level of an individual scattering. 
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